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We investigate the electronic states of the C0O2 plane in the layered cobalt oxides Nao.sCoCh 
using the realistic 11 band d-p model on a two-dimensional triangular lattice. Effects of the 
Coulomb interaction on a Co site: the intra- and inter-orbital direct terms U and U' , the 
exchange coupling J and the pair-transfer J', are treated within the Hartree-Fock approxima- 
tion. It is found that the metallic antiferromagnetism takes place below T c i and, in addition, 
the orbital order accompanied by the small charge order takes place below T C 2(< T c i) where 
the system becomes insulator. The obtained results are consistent with the successive phase 
transitions observed in Nao.sCo02. 

KEYWORDS: cobalt oxides, antiferromagnetism, charge order, orbital order, metal-insulator transition 
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The discovery of superconductivity in the layered 
cobalt oxide Na^OoC^-yEI^O 1 ) has stimulated further in- 
terest in the electronic states of the mother compound 
Na^OoC^. The specific features of the system are the ge- 
ometrical fluctuation of the C0O2 plane which consists 
of a two dimensional triangular lattice of Co atoms, and 
the orbital degeneracy of three t^g bands for Co 3d elec- 
trons. The electron filling in the C0O2 plane is controlled 
by changing the Na content x: the electron filling of the 
three ti g bands is given by 5 + x. For x < 0.6, Na^CoC^ 
is a normal paramagnetic metal except for x = 0.5, 
while, for 0.6 < x < 0.75, an anomalous behavior is 
observed: the magnetic susceptibility \ IS Curie-Weiss 
like although the resistivity p is metallic, 2 - ) the electronic 
specific heat coefficient 7 is large and increases with x. 3 ^ 
For x > 0.7, Na x Co02 exhibits large thermopower and is 
a good prospect for thermoelectric applications. 4 ~ 6 ' For 
x > 0.75, a weak magnetic order is observed below 20 
K, 7 ) where the ferromagnetic ordered C0O2 layers cou- 
ple antiferromagnetically with each other. 8 ) 

Furthermore, Nan.5Co02 exhibits remarkable succes- 
sive phase transitions at T c \ ~ 87 K and T C 2 ~ 53 
K: 2 ) the itinerant antiferromagnetism is realized below 
T c i and the system becomes insulator below T C 2- The 
magnetic susceptibility x shows kinks at T c \ and T c i. 
The resistivity p exhibits only a tiny anomaly at T c ±, 
while it rapidly increases below T C 2- From the NMR and 
the neutron measurements, Yokoi et al. 3 ^ have proposed 
the magnetic structure of Nao.sCoO^ below T cl , where 
chains of Co 3 ' 5 " 1 " 5 with larger staggered moments within 
the C0O2 (ab) plane and chains of Co 3 - 5- " 5 with smaller 
moments along the c axis align alternatively on the ab 
plane. The ordering of the c axis oriented moments on the 
Co 3 ' 5 " 15 chains has not been distinguished whether fer- 
romagnetic or antiferromagnetic, 3 ) and has not been de- 
tected by polarized neutron measurements. 9 ^ The charge 
disproportionation of the Co sites into chains of Co 3 5+l5 
and Co 3 5- " 5 is closely related to the ordered pattern of 
Na ions which form onedimensional zigzag chains below 
room temperature. 10-11 ' 



From the zero-field NMR spectra which exhibit split- 
tings below 7^2, 3 -* Ning et ai. 12 ^ have proposed the charge 
ordering pattern below T C 2, where the Co valence alter- 
nates on the Co 3 ' 5- * 5 chains while it is constant on the 
Co 3.5+«5 chairis . xhe 

zero-field NMR measurements have 
also been carried out for K0.5C0O2, 13 ' 14 ' which exhibits 
similar successive transitions at T cl ~ 60 K and T c2 ~ 20 
K, but the splitting below T c2 has not been observed in 
contrast to Nao.5Co02- Then, the origin of the metal- 
insulator transition at T C 2 still remains controversial. 

Many theoretical works have been done on the elec- 
tronic states of Narj.5Co02 for the last few years. 15 ~ 19 ) 
Some of the proposed magnetic structures are consis- 
tent with the experimental prediction below T c i, 3,9 ) how- 
ever, none of the proposed charge ordering patterns 
are consistent with the experimental prediction below 
T C 2- 12 ^ In our previous works, 20, 21 - ) the electronic states 
of Nao.5Co02 have been investigated using the 11 band 
d-p model on a two-dimensional triangular lattice, where 
the tight-binding parameters are determined so as to fit 
the LDA band structure 22 ) (see Fig. 1). Effects of the 
Coulomb interaction at a Co site are treated within the 
Hartree-Fock approximation, and the effect of the onedi- 
mensional Na order is considered by taking into account 
of an effective potential on the C0O2 plane. It was found 



TB with hole pockets (a) 
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Fig. 1. The band structure (a) and the Fermi surface (b) from the 
d-p tight-binding model with hole pockets (solid lines), without 
hole pockets (dashed lines) and from the LDA 22 ' (open circles), 
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that the Na order enhances the Fermi surface nesting, 
which is suppressed due to the frustration effect in the 
case without the Na order, and results in the metallic 
antiferromagnetism, 20 ) whose magnetic structure is con- 
sistent with the experimental prediction. 3,9 ) 

In addiction, we have considered the effect of the 
Coulomb interaction between the nearest-neighbor Co 
sites V and found that a coexistence state of the an- 
tiferromagnetic, charge and orbital order takes place for 
V > V c , where the system becomes insulator. 21 ) The 
obtained charge ordering pattern is consistent with the 
experimental prediction by Ning et aZ., 12 ) but the charge 
differentiation due to the charge order seems to be too 
large as compared to the splittings of the zero-field NMR 
spectra, which are very small for Nao.5Co02 and are ab- 
sent for K0.5C0O2. 13 ' 14 ) The overestimation of the charge 
differentiation is considered to be due to the effect of the 
intersite interaction V which is known to enhance the 
charge order. Therefore, in the present work, we discuss 
the spin, orbital and charge ordered states exclusively 
due to the effects of the on-site interaction in the ab- 
sence of the intersite interaction V. 

To investigate the electronic states of the C0O2 plane 
in the layered cobalt oxides NazCoC^, we employ the 
two dimensional triangular lattice d-p model 20,21,23,24 ) 
which includes 11 orbitals: d 3z 2_ r 2, d x 2_ y 2, d xy , d yzi d zx 
of Co and p Xx , pi v , pu (p 2x , p 2y , Pi z ) of O in the upper 
(lower) side of the Co plane. The model is given by the 
following Hamiltonian: 
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where dj^ (d„ m(T ) is a creation operator for a Co 
3d electron with wave vector k = {k x ,k y ) (site n = 
(n x ,n y )), orbital m(= 3z 2 — r 2 , x 2 — y 2 ,xy, yz, zx) and 
spin <r(=T) I), and p^ji a is a creation operator for a oxy- 
gen 2p electron with wave vector k, site j(— 1, 2), orbital 
/(= x,y,z) and spin a, respectively; 
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E 



t 



pd 



and t 



,1,1 



~ dnmo-dnma 

\ ma . The transfer integrals 

/, which are written by the Slater- 



Koster parameters, together with the atomic energies e p 
and e^j are determined so as to fit the tight-binding en- 
ergy bands to the LDA bands for Nao^CoC^ 22 -* as shown 
in Figs. 1(a) and (b). As the six hole pockets near the K 
points predicted by the LDA 22 ' have not been observed 
by ARPES experiments, 25 - 1 we employ the energy bands 
without the hole pockets where the tight-binding param- 
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Fig. 2. The schematic diagram of the antiferromagnetic ordered 
state (a) and the coexistence state of the antiferromagnetic, 
charge and orbital order (b). The filled (open) circles represent 
Co atoms on (out of) Na ordered lines and the upward (down- 
ward) triangles represent O atoms in the upper (lower) layers of 
the Co layer. The parallelogram is the extended unit cell. 



eters relating to the trigonal distortion are slightly mod- 
ified (see Figs. 1(a) and (b)). 

In the Hamiltonian eq. (1), we consider the effects of 
the Coulomb interaction at a Co site: the intra- and 
inter-orbital direct terms U and U' , the exchange cou- 
pling J and the pair-transfer J'. Here and hereafter, we 
assume the rotational symmetry yielding the relations: 
U' = U - 2J and J = J', and we set J = [7/10 for 
simplicity. We also consider the effect of the Na order 
at x — 0.5, where Na ions form onedimensional chains 
below room temperature, by taking into account an ef- 
fective onedimensional potential on the C0O2 plane 20,21 ) 



for odd n y 
(on Na ordered lines) 

for even n v 
(out of Na ordered lines) 



(2) 



with the effective potential Ae^ due to the Na order. In 
this paper, we set Aed = 1.0 eV. Due to the effect of 
Aed, the band structure becomes quasi one dimensional 
and the Fermi surface nesting is enhanced resulting in 
the antiferromagnetic order. 20,21 ) 

Now we discuss the possible ordered states of the 
Hamiltonian eq. (1) within the Hartree-Fock approxima- 
tion, where we assume that the order parameters are 
diagonal with respect to the orbital m and the spin a. 
Figure 2(a) shows the obtained ordered pattern of the 
antiferromagnetic state (AFM). Chains of Co 3 5+<5 and 
q q 3.5-a a ijg n alternately within the C0O2 plane due 
to the effect of the onedimensional Na order. 10 ) Fig- 
ure 2(b) shows an order pattern of a coexistence state 
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Fig. 3. The sublattice magnetization (n-f — n.jj/2, the orbital or- 



of states at the Fermi level Dp as functions of U for Ae^ 
eV at T = 20 K. 



1.0 



of antiferromagnetic, charge and orbital ordered state 
(AFM+CO+OO), where the Co sites out of the Na or- 
dered line with the valency of 3.5+<5 have the antifer- 
romagnetic moment, 3,9 ' while, the Co sites on the Na 
ordered line with the valency of 3.5-5 show a charge or- 
der of Co 3 - 5+s+r < and Co 3 - 5+<5 ~ 1 ' as proposed by Ning et 
al. 12 ) We note that the charge order in X-direction is dif- 
ferent from the charge order of Co 3 5+6 and Co 3 ' 5- * 5 in 
Y-direction due to the effect of the onedimensional Na 
order. 

Figure 3 shows the sublattice magnetization (nj — 
nj,)/2, the orbital order n 



sx , the charge order 77, 
the density of states at the Fermi level Dp and the re- 
sistivity p as functions of U at T — 20 K. Here we esti- 
mate the resistivity by using the Boltzmann equation as 
p = 1/ (a xx + <7 yy ) with 



— y 
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df 
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(p = x,y). (3) 



When U is larger than a critical value U c ~ 1.45 eV, we 
observe the coexistence state AFM+CO+OO schemati- 
cally shown in Fig. 2 (b). When U increases for U > U c , 
both the sublattice magnetization (nj — «j.)/2 and the or- 
bital order n yz — n zx increase, while the density of states 
Dp decreases and finally becomes zero for U > 1.7 eV. 
The resistivity p rapidly increases for U > U c and be- 
comes almost infinity for U > 1.7 eV. Due to the struc- 
ture of the orbital order (see Fig. 2 (b)), the charge order 
77 is induced by the orbital order although the amount of 
77 is vary small or zero as shown in Fig. 3. 

We plot the temperature dependence of the order pa- 
rameters together with the resistivity for several values 
of U in Figs. 4 (a)-(d). We observe the successive phase 
transitions at T cl and T c2 which increase with U . Below 
T cl , the sublattice magnetization (rif — n.jJ/2 becomes 
finite and the AFM is realized as shown in Fig. 4 (a). 
Below T C 2(< Tci) j the orbital order n. 



zx and the 

induced tiny charge order 77 become finite and the coex- 
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Fig. 4. The temperature dependence of the sublattice magnetiza- 
tion (n,j — ri|)/2 (a), the orbital order n yz — n zx (b), the charge 
order r\ (c) and the resistivity p (d) for Ae^ = 1.0 eV at U = 1.5 
eV (solid lines), U = 1.6 eV (dotted lines) and U = 1.7 eV 
(dashed lines). 



istence state of AFM+CO+OO is realized as shown in 
Figs. 4 (b) and (c). The sublattice magnetization is en- 
hanced due to the effect of the orbital order below T C 2 
as shown in Fig. 4 (a). In the case with large U = 1.7 
eV, the resistivity p shows a small anomaly at T cl , while 
it rapidly increases with decreasing T resulting the insu- 
lating grand state as shown in Fig. 4 (d). On the other 
hand, in the cases with small U — 1.5 and 1.6 eV, p 
shows kinks at both T c \ and T C 2 resulting the metallic 
grand state. We note that, for larger value of U > 1.8 
eV, T c i merges with T C 2 and the direct transition from 
the paramagnetic state (PM) to the coexistence state of 
AFM+CO+OO is observed. 

Finally, we show the U-T phase diagram includ- 
ing the PM, the AFM and the coexistence state of 
the AFM+CO+OO for Ae d = 1.0 cV in Fig. 5. 
The transition between the PM and the AFM takes 
place at T — T c \ 1 while the transition between the 
AFM (or the PM) and the AFM+CO+OO takes place 
at T = T c2 . For U > 1.49 eV, the system shows 
the successive phase transitions at T c \ and T C 2 as 
shown in Fig. 4. Remarkably, we also observe the reen- 
trant phase transitions of PM-AFM-PM for 1.45 < 
U < 1.47 eV, and the reentrant phase transitions of 
PM-AFM-PM-AFM+CO+OO for 1.47 < U < 1.49. 
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Fig. 5. U-T phase diagram for Ae,j = 1.0 eV. 



Such reentrant transitions are due to the schematic struc- 
ture of the density of states which have a minimum 
around the Fermi level under the one-dimensional po- 
tential Ae^. It is interesting to observe such reentrant 
transitions in experiments under high pressure. 

In summary, we investigated the electronic states of 
the C0O2 plane in the layered cobalt oxides Nao.sCoC^ 
by using the 11 band d-p model on the two-dimensional 
triangular lattice within the Hartree-Fock approxima- 
tion. What we have found are; (1) The system shows 
the successive phase transitions where the AFM state 
is realized below T c \ and the coexistence state of 
AFM+CO+OO is realized below T c2 (< T cl ). (2) The 
effect of the onedimensional potential due to the Na or- 
der is crucial for the successive phase transitions. (3) 
The resistivity shows a small anomaly at T c i, while it 
rapidly increases with decreasing T resulting the insulat- 
ing grand state for large U. (4) In the coexistence state 
of AFM+CO+OO, the charge order is induced by the or- 
bital order and is extremely small. The obtained results 
arc qualitatively consistent with the experimental results 
for the successive phase transitions in Nao.5Co02- 

In the present paper, we employed the energy bands 
without the hole pockets. In the case with the hole pock- 
ets, we also carried out the Hartree-Fock calculations. 
Due to the effect of the one dimensional potential, the 
resulting Fermi surface is almost the same as that for the 
case without the hole pockets. Therefore, the system with 
the hole pockets also shows the smiler successive phase 
transitions, although the AFM is slightly stabilized as 
compared to the case without the hole pockets. 

We have proposed a mechanism of the metal-insulator 
transition at T C 2, where the orbital order together with 
a small induced charge order takes place. The obtained 
charge ordering pattern is consistent with the experi- 
mentally proposed ordering pattern in Nao.5Co02, 12 - > al- 
though the orbital order, which is crucial in the present 
theory, was not considered there. There is another type 
of the coexistence state of the antiferromagnetism and 



the orbital order without the charge order, which seems 
to be consistent with the experiments in Ko.sCo0 2 13 ' 14 - ) 
where the charge order has not been observed below T C 2- 
The explicit calculation including this type of the coex- 
istence state is now under the way. 
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